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1.0 SUMMARY 


An improved integrated system of computer programs has been developed for the 
design and analysis of supersonic configurations. The system uses linearized 
theory methods for the calculation of surface pressures and supersonic area 
rule concepts in combination with linearized theory for calculation of 
aerodynamic force coefficients. Interactive graphics are optional at the 
user's request. 

The integrated system consists of an executive driver and eight basic computer 
programs including a geometry input module, which are used to build up the 
force coefficients of a selected configuration. 

The description of the design and analysis system is covered in four separate 
documents with the following subtitles: 

Part 1 - General Description and Theoretical Development (CR-3351) 

Part 2 - User's Manual (CR-3352) 

Part 3 - Computer Program Description (CR-3353) 

Part 4 - Test Cases (CR-3354). 

These documents supersede NASA Contractor Reports CR-2715, CR-2716, and 

CR-2717, which described an earlier version of the system. 

This document (Part 1) presents a general description of the system and 
describes the theoretical methods used. 






2.0 INTRODUCTION 


Over a period of years, NASA-Langley has developed a basic computerized series 
of supersonic design and analysis methods for aerodynamic configuration 
studies (reference 1). The methods are characterized by their reliability in 
use and input simplicity. 

The Boeing Company has extended this basic series of methods and combined them 
into an integrated system of computer programs. The extensions to the methods 
provide several new features: 

• Addition of a near-field (thickness pressure) wave drag program, to 
complement the existing supersonic area rule program. 

t Improved modeling of fuselage in lifting surface design and analysis 
programs. 

• Addition of configuration-dependent loadings in the wing design 
program, so that the wing design is performed in the presence of 
fuselage and nacelle effects. 

• Addition of pressure limiting terms in the lifting pressure programs, 
to constrain the linear theory solution. 

• Optional CRT displays of selected program input and output data and 
provisions for limited user editing and intervention. 

A plot module is included in the system to produce configuration drawings, and 
a common geometry module is used to permit a single geometry input for all 
programs. A wing pressure module permits summaries of wing pressures at 
desired conditions. 

The basis of the system is supersonic linearized theory, modified in two 
respects: 

• The "Whitham" correction to disturbance positioning is used in the 
propagation of body pressure fields. 

• The wing lifting pressure modules contain an optional limiting 
pressure feature to control the permissible level of upper surface 
pressure coefficient. In addition, the wing design module contains a 
further constraint feature to limit the upper surface streamwise 
pressure gradient. 

Superposition is used to build up the theoretical force coefficients of a 
selected configuration. 

The goals of the integrated system have been to develop an easily used 
supersonic design and analysis capability, with recognition of the need for 
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constraints on linear theory methods to provide physical realism, and with 
inclusion of interactive display for increased design control over 
optimization cycles. 

The description of the design and analysis system is broken into four 
documents: 

Part 1 - General Description and Theoretical Development (CR-3351) 

Part 2 - User's Manual (CR-3352) 

Part 3 - Computer Program Description ( CR- 3353) 

Part 4 - Test Cases (CR-3354). 

This report (Part 1) presents a general description of the system and the 
theoretical methods used. 

These four documents supersede NASA contractor reports CR-2715, CR-2716, and 
CR-2717, which described an earlier version of the system. 
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3.0 SYMBOLS 




A 


A (x) 

B 


C 

C 




mo 


c p> Tim 


element fraction 
load strength factor 
cross-sectional area at x 
trail ing-edge element fraction 

pitching moment reference length 
wing tip element fraction 
chord 

section axial force coefficient 

section drag coefficient 

drag coefficient, D/qS 

minimum drag coefficient 

skin friction coefficient 

section lift coefficient 

lift coefficient, lift/qS 

lift coefficient per degree angle of attack 

pitching moment coefficient, pitching moment/qSc 

pitching moment coefficient at zero lift 

section normal force coefficient 

local pressure coefficient, (p - p^ J/q^ 

limiting pressure coefficient used in definition of 
attainable thrust 
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ACp 

(ac p v^t 0 

c t 

★ 

c t 

C T 

D 

F(y) 

h(Z) 

k 

K i> k 2 , k 3 
k t 

L 

M 


P 

P< 

q 

R 

R 

r 

S 


lifting pressure coefficient (lower surface pressure 
coefficient minus upper surface pressure coefficient) 

leading-edge singularity parameter 

limiting value of leading edge singularity parameter at 
x* * 0 

section thrust coefficient 

attainable section thrust coefficient 
theoretical wing thrust coefficient 
drag force 
Whitham function 

decay function for F(y) calculation 

Mach number parameter 

constants used in curve-fitting 

attainable thrust factor, c t */ c t 

grid element variable in x (lengthwise) direction 

Mach number 

grid element variable in spanwise direction 
pressure 

freestream static pressure 

2 

dynamic pressure, 0.7 p w M 
body radius 

influence function 

radial distance to point in flow field 

reference area or cross-sectional area (figure 4.3-3) 
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T 


t 

U 

u 

v 

V 

W(Z) 

w 

X 

X I 

y 

z 

Z 

a 

e 

B 

£ 

X 

n 

p 

4 > 

9 

A 


temperature 

variable of integration (figure 4.3-4) 
freestream velocity 
x perturbation velocity 
y perturbation velocity 
freestream velocity vector 
body pressure decay function 
upwash, w perturbation velocity 
lengthwise variable 

distance behind local wing leading edge 

span station or distance variable 

vertical variable 

decay function parameter 

angle of attack 

local fuselage camber angle 

Mach number parameter = yM^ - 1 

x variable of integration 

surface shape (equation 13) or Lagrange multiplier 
(figure 4.4-4) 

spanwise variable of integration or semispan fraction 
density 

velocity potential 
radial angle variable 
leading-edge sweepback angle 
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T 

M 


Subscripts: 

oo 


i 

U j 

t 

Superscripts: 


ii 


* 


angle at a cambered-wing leading edge between tangent to 
surface and wing reference plane 

region of influence 

dynamic viscosity coefficient 


freestream conditions 

incompressible condition (equation (5)) 

loading numbers (equation (31)) 

tilted lengthwise variable (figure 4.3-5) 



reference temperature condition (equation (2)) or. first 
derivative with respect to x 

second derivative with respect to x 

field point grid element 
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4.0 DISCUSSION 


The integrated supersonic design and analysis system is shown schematical ly in 
figure 4.0-1. It consists of an executive "driver", and eight basic computer 
programs. The individual programs, or modules, provide data for configuration 
design or analysis as follows: 

t Skin friction is computed using turbulent flat plate theory. 

• Wave drag is calculated from either far-field (supersonic area rule) 

or near-field (surface pressure integration) methods. The far-field 
method is used for wave drag coefficient calculations and for 

fuselage optimization according to area rule concepts. The 
near -field method is used primarily as an analysis tool, where 
detailed pressure distributions are of interest. 

t Lifting pressures, drag-due-to-1 ift, pitching moment and trim drag 

are computed from the lift analysis program., which, breaks arbitrary 
wing/fuselage/ canard/nacelles/horizontal tail configurations into a 
mosaic of "Mach-box" rectilinear elements which are employed in 
linear theory solutions. A complementary wing design and 

optimization program solves for the wing shape required to support an 
optimized pressure distribution at a specified flight condition. 

• The plot module draws configuration pictures according to size and 
view inputs. 

• The geometry module handles configuration geometry for the system. 

• The wing pressure module summarizes and tabulates wing surface 
pressure data for input conditions. 

t Interactive graphics, for data display or editing, are optional in 
the system at user’s request. 

Operating in the analysis mode, the force coefficients of a selected 
configuration are built up through superposition as shown in figure 4.0-2. 

The entire design and analysis system is a single overlaid program. The 
executive level of the system controls module execution by means of special 
identification cards in the input data. Transfer of data in the system 
between modules is handled by disk storage and common blocks. 

All configuration geometry is read (or updated) in the geometry module. 
Geometry inputs are patterned after those of the NASA- LRC configuration 
plotting program. All "paneling" of the configuration for theoretical 
analyses is accomplished within the programs, and the user prepares only 
"drawing-type" geometry. 
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Zero Lift Wave Drag (Wave and vortex) 


FIGURE 4.0-1.— integrated supersonic design and analysis system 







•n di ^? a V a \ a ° n i ys , tem structure af id input and output formats 
in Parts 2, 3, and 4. The theoretical methods used in the 
in the remainder of this document. 


are presented 
modules are 


4.1 Skin Friction Drag 

m^hnH 1n n/ ri ' C V* 0n dra 9 of arbitrary configuration is calculated from the T' 
method of reference 2. Smooth flat-plate, adiabatic-wall, and turbulent 

leading J^ er f COnd h tl0n |* are assumed » with transition assumed to occur at the 
leading edge of each configuration component. 

ill* r“^rence nd 2 eXPe {Vlr en ™ verification of the T' method are given in detail 

tonendifr Of rL Jr® eSSent \ a1 elements of the method are presented in 
Appendix C of the reference, and are summarized here for completeness. 

coefficient 0 ^ 1S f baSed ° n l he calculation of a compressible skin friction 
coefficient, Cp, from a reference skin friction coefficient C'r for a 

temperature 31 "? "“"mf M h “ ’■ , Reyno ! ds number, R„ , and adiabat’ic wall 
conditions).’ W ' subscnpt - is used throughout to denote free stream 

^lationlh,-nf«c, atUre rat l?’ V T - ’ is calculated from one-dimensional 
relationships assuming a wall recovery factor of .89. 


W 2 

■ 1 + 0.178 Mi 


(i) 


relationship f ° W * S ° mmer a " d Sh ° rt 1n reference 2 obtained the kej 
■f 1 = 1 + 0.035 K„ 2 + 0.1*5 Q*-- l) 

( 2 ) 


The Reynolds number relationship is 

1 


R' 

R_ 


©(© 


(3) 


where the viscosity ratio is given by the Sutherland equation 


M 1 _ / T' V ( T «o + \> 

P- " \ Too/ '' TATIS'/ 


(4) 


with T in degrees Rankine. 

The incompressible skin friction coefficient is given by the Karman-Schoenherr 
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equation 


0.242 


1Og 10 ( C F. R i) 


(5) 



V c 

F i 


which gives the corresponding relationship based on the T' analogy: 


0.242 

vy 


log 


10 


( 


0 


( 6 ) 


This equation is solved iteratively for Cp'. The desired compressible skin 
friction coefficient is obtained from Cp' by 


C 


F 




(7) 


Most of the skin friction program is involved with computing wetted areas and 
reference lengths for the configuration. Components which may exhibit large 
variations in reference length (such as the wing, tail, canard etc. ) are 
broken into strips in order to calculate Cp accurately. 



4.2 Far-Field Wave Drag 

The far-field wave drag program computes the zero lift wave drag of an 
arbitrary configuration by means of the supersonic area rule. A description 
of the program has been given in reference 3, and is summarized here. The 

method assumes not only an integrated total lift of zero, but local lift 
everywhere zero, and thus includes none of the wave drag contribution 
associated with the generation of lift. 

The supersonic area rule is a generalization of the transonic area rule, which 
states that the transonic wave drag of an aircraft is essentially the same as 
the wave drag of an equivalent body of revolution having the same 
cross-sectional area distribution. In the supersonic area rule procedure, 

several equivalent bodies of revolution are produced by passing a series of 
parallel cutting planes through the configuration, as shown below. 
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The cutting planes are inclined with respect to the aircraft axis at the Mach 
angle p , and a single equivalent body is produced for the series of cutting 
planes at a constant azimuthal angle, 0. The area of the equivalent body at 
each station is the projection of the area intercepted by the cutting plane 
onto a plane normal to the aircraft axis. 

The wave drag of each equivalent body is calculated from the Von KaVman 
slender-body formula (reference 4), which gives the drag as a function of the 
free stream conditions and equivalent-body area distribution. 

pu 2 ft (8) 

D(e) = ““mT J J A " (x 1 } A " 106 I X 1 “ X 2^ ^1 ^2 


where xi and X2 are lengthwise variables of integration and 
A" is trie second derivative of the body area distribution. 

The wave drag of the aircraft at a given Mach number is calculated from the 
integrated average of the equivalent-body wave drags. 

1 ( ^ (9) 

D * 2T / D(e) d0 

J o 

A useful feature of the supersonic area rule occurs in the optimization of a 
fuselage area distribution to minimize the wave drag of a wing-fuselage 
combination. According to Sheppard (reference 5), the wave drag of a 
wing-fuselage can be written as: 


-HI) *(i) 


wing eq. body 
+ body 


-(i) 


where the wing equivalent body is a body of revolution obtained by averaging 
from 0 = 0 to 2tt all the projected areas intercepted by cutting planes 
passing through the wing for each X station. 


Considering a wing-fuselage, then, and assuming that the wing geometry is 
fixed, both the first and third terms of equation (10) are fixed. This leaves 
optimization of the second term only, so that to design the minimum drag 
configuration, the fuselage must be contoured to produce a wing-equivalent- 
body-plus-fuselage having the shape of a minimum drag body of revolution. 


The wave drag program uses the Eminton-Lord fairing through a discrete set of 
points, which defines the shape of a minimum drag body through the points. 
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Using this fairing, and identifying a few "control-points" on the fuselage 
area definition, an optimized fuselage area definition may be obtained as the 
difference between the wing equivalent body and the combined wing-equivalent 
-body-plus-fusel age. 

The wave drag program has mechanized this solution, and also includes the 
effects of nacelles, tail, etc., which are included similarly. The shape of 
the minimum drag fuselage area distribution, and the drag of the complete 
configuration including the minimum drag fuselage, may be calculated for each 
configuration input. 

Because of the configuration generality that can be handled, the far-field 
program is the primary source of zero-lift wave drag coefficients in the 
design and analysis system. A complementary near-field program, more 

restricted in scope but which computes detailed thickness pressures, is 
described in the following section. 


4.3 Near-Field Wave Drag 


Isolated Component Pressures 

The near-field wave drag program computes zero-lift thickness pressure 
distributions for wing-body-nacelle-empennage configurations in supersonic 
flow. The pressure distributions are integrated over the cross-sectional 
areas of the configuration to obtain the resultant drag force. Three basic 
calculations are performed to obtain the required pressure fields: 

• Thickness pressure distribution for a wing or tail of arbitrary 
section and planform. 

• Thickness pressure distribution on surface of fuselage or nacelles. 

• "Whitham" near-field calculations to define pressure distributions 
propagating from fuselage or nacelles. 

Superposition is used to calculate the interference drag terms associated with 
the pressure field of a given component acting on the surfaces of the other 
components. 

The surface pressure distributions calculated include some (from the nacelles) 
which contribute lift due to a non-symmetrical distribution of volume, so that 
the drag calculations which include these loadings contain wave drag due to 
lift contributions. However, by selection of program options and subtraction 
of certain drag contributions, a zero lift wave drag comparable to the 
far-field wave drag value may be obtained, as described on page 34. 


15 



Wing or empennage thickness pressures . - The surface pressure coefficient on 
the upper (or lower) surface of a flat mean-line wing of symmetrical surface 
shape is obtained by first calculating the corresponding velocity potential, 
differentiating with respect to x (to get u), and then computing the pressure 
coefficient from the linear theory approximation, C p = -2U/U*, . 

The velocity potential computation, from reference 13, is 


<t> (x,y) 


1 if X dn dg 
¥ JJ^x-c) 2 - B^y-n) 2 

T 


(ID 


where $ (x,y) 
X 

B 

Z 

n 

T 


velocity potential at a defined wing field point (x,y) 
surface slope (dz/dx) of wing section at a wing 
in tegrat ion point 
VW~TJ 

x variable of integration 
y variable of integration 

subscript denoting interval of integration (surface of 
the wing planform within the Mach forecone from x,y) 


The wing thickness pressure coefficient is 


C p (x,y) = P— Pqp = -2 M = -2 3$ (x,y) 

Poo ^oo 3x 


( 12 ) 


where p = local pressure at x,y 

p = free stream static pressure 

u 00 = x perturbation velocity 

^oo >Poo = free stream velocity and dynamic pressure 

Integration of the velocity potential equation is performed by representing 
the wing as a grid of rectangular elements, and substituting a numerical 
summation for the integration. The wing element "Mach-box" system exactly 
follows the arrangement described in Section 4.4 (for the lifting pressure 
case) and is summarized in figure 4.3-1. The grid elements, identified by L 
and N, are defined such that L is numerically equal to x and N is numerically 
equal to By, where x and By take on only integer values. Partial grid 
elements along the wing leading and trailing edges are used to permit a closer 
approximation to the actual wing planform. The grid system of figure 4.3-1 is 
rather coarse for illustrative purposes; in actual usage, many more grid 
elements are employed. 

With respect to a specified field point x,y, the upstream region of 
influence! is approximated by the shaded grid elements in figure 4.3-1. Each 
of the shaded elements has associated with it an influence factor, R, which 
relates the effect of the element and its surface slope to the total velocity 
potential at the field point. 
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The factor R is obtained by an approximate solution to Equation (11) as 
follows: v ' 


Let y-n = u so that du = -dr\ 




(x,y) -4 




du d £ 


_ \ 2 2 2 

x-C) - B u 


(13) 


Substituting an average value of (x -£) - (x -£) and integrating with respect 
u 0 u , 


Restoring y - n = u 


4> (x,y) = — 

IT 


<t> (x,y) 


= -i-[ 

ttB L 



Bu 

(x^) 


X sin 


B(y-n) 

#rr 



(14) 


(15) 


For a given grid element in the Mach forecone, the integration interval rii 
to r )2 gives 1 

0 ( Y - r ) ) = N* - N + 0.5 

0{y-n 2 ) = N* - N -0.5 

Also, = L* - L + .05 

L* and N* define the field point grid element and L and N define an element in 
the forecone. The numerical summation then becomes 




L,N 



N* - N - 0.5 
L* - L + 0.5 


sin 


-1 


N» - N + 0.5 ~| 
L* - L + 0.5 J 


(16) 


with R being the quantity between the brackets for the element (L,N). 

The character of the R function is such that the sum of R for all elements 
across the Mach forecone at a constant L value is equal to -tt . In 
considering the case of an infinite rectangular wing of constant thickness 
shape (wedge section), it is seen that the -tt summation must be the case, 
since the $ function for successive L* rows is increased by 



L=*const . 
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so that 


X 



A *. . 

Ax 



and 



which is the established linear theory result. Values of R within the Mach 
forecone from a typical field element are presented in figure 4.3-2. 


The computer program for the wing thickness pressure calculations sets up the 
winq qrid system, defines the surface slope X for each grid element performs 
the summation for <f> at selected field elements, and differentiates 9 
numerically to obtain C p . Wave drag is obtained by summing the pressure 
coefficients times the corresponding surface slopes of the wing. 


In the wing surface slope definition, slopes of the elements occupying space 
covered by the fuselage (if there is one) are set equal to zero to eliminate 
their drag contribution. Partial elements are used at the wing-body 
intersection to improve the geometry definition. 

As was found to be the case in lifting surface programs involving the grid 
system used in the wing calculations, some oscillation in computed pressure 
coefficient values occurs. This has been suppressed by a smoothing of the 
calculated pressure coefficient values, using a 5 term equation of the form; 


.333 A 


L*-2 


C =■ 
P 


p L*-2 


+ .667 ^ 


P L*-1 


+ A. 


L* p 


+ .667 C 


+ .333 C 


( 17 ) 


L* 


L # +l 


\ L*+2 


333 + .667 A l# _ 1 + A l# + 1 


.0 


where A = factor defining size of associated element (1 if whole element, 
fraction if not). 



-uselage and nacelle thickness pressures . - The pressure distribution on the 
surface of fuselage or nacelles is obtained from a method developed at the 
Boeing Company by R. M. Kulfan, based on the Lighthill theory of reference 6. 
The method is applicable to bodies having either smooth area distributions or 
bodies with slope discontinuities. As used in the near-field program, smooth 
area distributions are assumed except (if required) at the nose or aft end of 
the body. Open-nose bodies, such as nacelles, are permissible. 


The solution technique is summarized in figure 
velocity is calculated which is a function of 


4.3-3. An axial perturbation 
the body cross-sectional area 
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growth (and radius distribution) and a decay function which relates area 
growth to its effect on a given field point. 

The axial perturbation velocity, u, is given by 


where 


x 



X = 

W(Z> = 

K = 
Z = 


: 

<*) 


body field station 
decay function 
x variable of integration 
position function = x ~ E 


body radius at £ 
first derivative 
area (S) at £ 



of body cross-sectional 


(18) 


The pressure coefficient is calculated from one-dimensional flow relationships 
(rather than the simpler linear theory approximation), as shown in figure 
4.3-3. The decay function, W(Z), is tabulated on page 23. 


Fuselage and nacelle flow-field pressures . - The pressure field propagated by 
the fuselage or nacelles is calculated using the Whitham solution of reference 
7. This solution, which is converted to a method amenable to the digital 
computer in reference 8, calculates shock wave positions and strengths through 
a modification of linear theory results. 


The Whitham solution begins with the calculation of a function, 
dependent on body geometry: 


F(y) 




1 

2 


h <*) a s- (t) 


F(y) , which is 


(19) 


where 

y = body field station 

h (Z) = decay function (similar to body pressure 
calculation) 

t = y variable of integration 
Z = t position function 

s'(t) = first derivative of body cross-sectional 
area (S) at t 

^end = Y at end integration (see figure 4.3-4). 

The solution for F(y) is very much like that of the body thickness pressure 
solution, with the exception that the integration is carried out to the point 
(y en d) at which the aft-running Mach line from the body centerline at y 
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0 R d> 


S = cross-sectional area 


where: 

dS may include dS smoot h ar| d dS discontinuous 


u(x) = Au (x) 

smooth 
AU '< x) disc = 



d£ 


Where integration 
excludes ds' disc 





t 

\ 

3.5 ) 

PRESSURE 

CpM 

' 2 

1.0- 0.2 ( 

(l.° + R'(x|^(l -0 + u |x |)^ -’.o) 


COEFFICIENT 


0.7 Mj I 

L \ 

/ J 

) 


FIGURE 4.3-3.-BODY SURFACE PRESSURES (LIGHTHILL THEORY) 



DECAY FUNCTION 




emerges from the body contour. The F(y) calculation is summarized in figure 
4.3- 4, and the decay function h(Z) is tabulated on page 27. 

Development of the near-field pressure signature at a radial distance (r) from 
the body centerline is summarized in figure 4.3-5. The F(y) function is 
tilted, as shown, by displacing the F(y) function according to 

Ay = k/r* F(y) (20) 


where , 

k _ 1.6973 M 4 

' 6 1 - 5 

This process results in a double-valued pressure signature where shock waves 
are present, a situation which is corrected by an area-balancing technique. 
The area-balancing (shaded lobes are equal in area) defines shock wave 
locations and strengths. The resultant tilted/area-balanced signature is then 
converted to pressure coefficient by 

c - P ~ P “ = a F(y) (21) 

P V 28r 

The position of the pressure signature at radius (r) is then given by 


x = y t + Br (22) 


where x is the longitudinal distance from the body nose. The "tilting" thus 
produces a correction in pressure signature positioning, relative to simple 
translation along a Mach line, which is the result of the remarkable theory of 
reference 7. 


Interference Calculations 


The previous sections summarized the basic pressure calculations performed by 
the near-field wave drag program. Integration of the wing and fuselage or 
nacelle pressure distributions over their respective surfaces gives the 
"isolated" wave drags of the components. Interference drags are obtained by 
superposition, i.e., the pressure field of each component is imposed upon the 
surfaces of the v other components to calculate the resultant interference 
forces. Interference terms between empennage surfaces and the rest of the 
configuration have been considered small, however, and are neglected. 

The computer solution allows for up to 3 pairs of nacelles located external to 
the wing-fuselage (or 2 pairs plus a single nacelle at y = 0). The nacelles 
may be either above or below the wing (or both). 
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DECAY FUNCTION 



PRESSURE P-Poo _ 2 F(v) 

COEFFICIENT r " 


where: 

r = radial distance from body 
F{y) function displaced, shock waves 
located as shown in fig. 4.3-5 



FIGURE 4.3-4.- WHITHAM F(y) FUNCTION (NEAR-FIELD PRESSURE BASIS ) 



k - 


1.6973 M 4 
j3 1,5 


, r = radial distance from body 
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Wing-fuselage interference . - Fuselage-on-wing interference is obtained by 
calculating the near-field pressure signature from the fuselage at selected 
spanwise stations, and imposing them upon the corresponding wing sections as a 
buoyancy field. The spanwise stations are the same as those used in the 
isolated wing thickness pressure calculations. 


Wing-on-fusel age interference is obtained from wing thickness pressures 
computed in the area occupied by the fuselage, after the wing surface slopes 
in the fuselage area have been set equal to zero. The thickness pressures 
thus calculated are "carry-over" pressures, and are imposed upon the fuselage 
surface slopes by transposing them aft along Mach lines to the fuselage. The 
fuselage area covered by the wing root is deleted from the interference term. 


Typical results from wing-body calculations are presented in figure 4.3-6. 

The wing data, both isolated and fuselage-on-wing interference, are converted 
to section drag form. 


Nacelle interference terms are calculated like the wing-body interference 
terms, with the provision that two alternative solutions may be specified for 
the nacelle-on-wing interference term. Available experimental data do not 
make it clear whether a "wrap" or "glance" solution, as shown in figure 4.3-7, 
is more correct. Since the nacelle-on-wing interference term is substantial' 
both solutions are available in the program (controlled by an input code). 

The nacelle interference terms are calculated as follows: 


N jc^ile-on-wlng term . . - The nacelle pressure field acting on the wing is 
obtained by calculating nacelle pressure signatures at the same spanwise 
stations used for the wing thickness pressures (plus extra stations 

immediately adjacent to the nacelle centerlines), then defining a composite 
signature by summing the contributions from all nacelles. The nacelle 
pressure coefficients are doubled to account for reflection from the winq 
surface. 

In the case of nacelles both above and below the wing, separate nacelle 
pressure fields for nacelles below, and then above, the wing are calculated 
and then sumned. 

A unique feature of the near-field approach lies in the solution choices 
available, since the pressure signatures generated by a nacelle to act on the 
wing surface may be terminated on encountering another nacelle or may be 
allowed to pass undiminished around (or through) it. In the first case, the 
'glance" solution, the nacelle pressure field is assumed to reflect from other 
nacelles in such a way that the reflected field exerts no further influence 
upon the wing surface. In the second case, the "wrap" solution, the nacelle 
pressure field is assumed to propagate around the other nacelles as if no 

obstruction were offered; i.e., as if they were transparent to the pressure 
field. 
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FIGURE 4.3-6. - WING-BOD Y SOL UTION , M ^ 2.6 



PRESSURES "GLANCE” AWAY FROM WING AT ADJACENT 
NACELLES 

The nacelle pressure field and accompanying shock waves "glance" 
away from the wing when encountering adjacent nacelles. In 
application, the nacelle generated pressure field is terminated on 
encountering another nacelle. 
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PRESSURES "WRAP" AROUND ADJACENT NACELLE 

The nacelle pressure fields and accompanying shock waves 
"wrap" around adjacent nacelles. In application, the nacelle 
generated pressure field is allowed to pass through another 
nacelle as if it were transparent. 


FIGURE 4.3-7.— NACELLE GENERA TED WING PRESSURE FIELD CONCEPTS 




V = 12.9 




© Experimental data 

— Theoretical solution (wrap or glance) 



FIGURE 4.3-8.-NACELLE PRESSURE FIELD DATA , M = 2.7 


A comparison of the theoretical nacelle pressure field with experimental data 
is^ presented in figures 4.3-8 and 4.3-9 at Mach numbers of 2.7 and 1.1. At 
M = 2.7, there is. no difference between wrap and glance solutions. At 
M = 1.1, the solution choice appears to favor the glance approach except on 
the outboard wing. However, the experimental pressure points are too sparse 
to permit firm conclusions. 


Nacel le-on-fuselage term . - The effect of the nacelles on the fuselage area 
distribution is built up by integrating each nacelle's pressure signature upon 
the fuselage area growth, with the provision that nacelles below the wing 
affect only the fuselage area below the wing in the region of the wing (and 
vice-versa for above-wing nacelles). In the fuselage region below the wing, 
or above the wing, the nacelle pressures are doubled to account for reflection 
(as is done for the nacelle pressures acting on the wing). 


Fuselage-on-nacelle term . - The fuselage pressure signature is imposed on 
each nacelle surface as a buoyancy force. 


Wing-on-nacelle term . - Wing thickness pressures at span stations where 

nacelles are located are transferred aft along Mach lines from the wing 
surface to the nacelle centerline to obtain the buoyancy field acting on each 
nacelle. 


Nacelle-on-nacel le term . - The interference term of other nacelles acting on 
a selected nacelle is calculated by building up the composite buoyancy field, 
and then imposing it upon the nacelle surface. In the case of nacelles on 
opposite sides of the wing ( i . e . , above and below wing), the pressure 
signatures are cut off where intercepted by the wing. 

The nacelle-generated pressure field is assumed to pass undiminished around or 
through other nacelles that may be present when a particular nacelle pair 
interference term. is being calculated much as in the "wrap" solution employed 
for nacelle on wing terms. The "wrap" or "glance" option is not provided in 
this case because of the generally negligible difference in results. 


Nacelle image effects . - If the nacelle is located next to the wing (or 
body), an "image" nacelle is used to create the reflected buoyancy field. The 
reflected field is cut off forward or aft of the wing if the nacelle is not 
entirely under the wing. 

The principal image effect is caused by the nacelle pressure field reflecting 
off the wing back onto itself. However, the reflected nacelle pressure field 
acting on other nacelles is also computed. This solution also utilizes an 
image nacelle representation, and the reflected signature is restricted to 
whatever part of the "real" signature encountered a reflecting surface. 
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Composite Thickness Pressure Signatures and Configuration Drag 


Typical theoretical thickness pressure distributions acting on the 
configuration components are presented in figure 4.3-10. The net effect of 
these pressure signatures on drag is obtained by summing the various inputs 
and integrating them over the -surface of the configuration. The corresponding 
drag terms are summarized in the data of figure 4.3-11. 

The near-field wave drag program has three calculation features that have no 
direct counterpart in the far-field wave drag program. These are: 

1) Wrap or glance solution for nacelle pressure field acting on wing. 
The wrap solution is nominally the same as the far-field solution. 

2) Nacelle image effects. Because of the transparency assumption of the 
far-field program pressure propagation, nacelle pressures do not 
reflect off adjacent surfaces. (Through addition and subtraction of 
separate calculations involving "image" nacelles, however, comparable 
drag data can be generated using the far-field program.) 

3) Above- or below-wing fuselage area separation in the nacelle-on- 
fuselage term. Directly comparable results are obtained only for 
mid-wing configurations. 

Using the "wrap" solution, a mid-wing arrangement, and subtracting the nacelle 
image drag terms, a direct comparison between near-field and far-field proqram 
drag calculations can be made. 


4.4 Drag-Due-to-Lift (Design and Analysis) 


The wing design and lift analysis modules are separate lifting surface methods 
which solve the direct or inverse problems of: 

• Design - to define the wing camber surface shape required to produce 
a selected lifting pressure distribution. The wing design program 
includes methods for defining an optimum (least drag) pressure 
distribution. 

• Lift, analysis - to define the lifting pressure distribution acting on 
a given wing camber surface, and calculate the associated force 
coefficients. 




FIGURE 4.3-10.- NEAR-F/ELD PRESSURE SIGNATURES, M = 1.1 "WRAP" SOLUTION 
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Originally, these programs were developed to utilize the "Mach-box" method of 
wing representation and evaluation of linear theory integral equations 
(references 9 through 11). In the Mach-box method, the winq is replaced by a 
grid of small rectilinear elements. Since many elements (in the thousands) 
can be used, a detailed description of complicated surface shapes, with 
associated computational accuracy, is possible. 

The wing design and lift analysis programs have been expanded several times to 
add features (e.g., reference 10) and have been reviewed in reference 12. The 
discussion of the aerodynamic theory of the Mach-box method in reference 12 is 
quite detailed, and is summarized here for completeness. The numerical method 
for the "direct" case of wing design and optimization is given first, followed 
by the inverse case of lift analysis on page 53. 


Design Case 


Camber surface for a given loading . - A typical wing planform described by a 
rectangular Cartesian coordinate system is illustrated in figure 4.4=1. — The y 
coordinate has been multiplied by the Mach number parameter 0 =yM^ - 1 for 
convenience in mathematical manipulations. 

For a wing of zero thickness lying essentially in the z=0 plane, linearized 
theory for supersonic flow defines the wing surface shape necessary to support 
a specified lift distribution by the integral equation 

(x - O AC U,n) dn 
P . . — 

( y - 7i) 2 1 /(x - Tf - e 2 (y - n) 2 


dZ 
C 

3x 


(x,y) = tj- AC (x,y) + 


■fc H. 


which is a slightly modified form of equation (77a) of reference 13. With 
respect to a specified field point x, y, the upstream region of influence, i, 
enclosed by the Mach forecone is shown by the shaded area. 

Equation (23) may be rewritten into an influence function form by introducing 
the factor R, such that 

(x,y) = =£-AC p (x,y) +-|j j=A5 


. 


R(x-e>y-n) 


AC U,n) 
P 


dBn 


(24) 


where R is defined as 

. x - K 

R(x-£,y-n) * ~~~~ • 0 ^ 

e 2 (y - n) 2 V (x - O - p (y - n) 


(25) 
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R is a function relating the local loading at points , to its influence in 
determining the downwash (or upwash) at the field point. 

For the numerical evaluation of this integral, the wing is replaced by the 
grid system of figure 4.4-2. The grid elements, identified by L and N, are 
defined such that L is equal to X and N is equal to By, where X and By take 
on only integer values. Partial grid elements along the wing leading and 
trailing edges are used to permit a closer approximation to the actual wing 
planform. The grid system of figure 4.4-2 is rather coarse for illustrative 
purposes. In actual usage, many grid elements are employed. 


In the grid system, the field point x, y is located at the aft midpoint of a 
field point element L*, N*, and the region of integration, x, is approximated 
by the shaded grid elements. 


The integration required by equation (24) is performed for each element within 
the Mach forecone, considering AC p constant over the element and using an 
average value of X- £ = L_*-L+0.5. The resulting expression, derived in 
reference 9, is the factor R which is the value of the influence function for 
an element L, N: 


R(L*-L,N*-N) 




» - L + Q.5) 2 - (N* - N - Q.5) 2 
(L* - L + 0.5)(N# - N - 0.5) 


(26) 


V 


(L» - L + 0.5) 2 - (N» - N + 0.5) 
(L* - L +-0.5)(N* - N + 0.5) 


2 


The integral equation (23) may then be replaced by the summation given below, 
where the summation includes all elements within the forecone and on the wing 
planform. The factors A, B, and C are element fractions for the wing leading 
edge, trailing edge, and wing tip, respectively, to allow for partial elements 
at those locations. 

3z (27) 

— -S-(L«,N«) = 4- AC (L»,N«) 

d X ** it 

B Y* R(L*-L,H*-N) A(L,N) B(L,N) C(L,N) AC (L,H) 

P 

T 


The character of the R function is such_that, for a given L*-L set of elements 
within the forecone, the sum of the R values is zero, the single negative 
value at N* - N = 0 balancing all the others (see figure 4.4-3). At L* = L, 
where there is only one element in the summation, the R value is zero. 
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The physical significance of this R variation is that (for positive lift), 

all elements directly ahead of the field point element contribute only 
downwash and all other elements contribute upwash. An element at the leading 
edge near the wing tip of a subsonic leading- edge wing, therefore, sees a 
concentrated upwash field. It is this upwash field that makes the subsonic 
leading edge twisted and cambered wing attractive from the standpoint of 
drag-due-to-1 ift, since a local element may be inclined forward to produce 
both lift and thrust. 

Conversely, an element located at the trailing edge and centerline of a 
diamond planform produces lift inefficiently, since the element sees an 
influencing region that contributes downwash. 

In computing the local surface slope required to obtain a specified lifting 
pressure coefficient, the Cp value is known everywhere on the wing from 
loading formulas. In the wing design computer program, the wing camber 
surface shape is found by integrating the surface slopes for all chordwise 
elements at selected span stations. Section values of drag ( Cp times slope), 
lift and pitching moment are then integrated spanwise to obtain wing force 
coefficients. 


Optimum combination of loadings . - Equation (27), which defines the surface 
slope to support a specified lifting distribution, requires a description of 
the design pressure distribution. This is obtained from Lagrange's method of 
undetermined multipliers, which provides a means of selecting an optimum 
combination of component loadings, yielding minimum drag subject to various 
imposed constraints on the allowable pressure distribution. 

A number of design point options are provided in the wing design program, 
involving the selection of loadings to be combined and the constraints 
desired. (The options are controlled by input codes, as described in the 
User's Manual, part 2). 

The usual design conditions specified are: 

• Drag-due-to-1 i ft of the wing be minimized at a given design lift 
coefficient, subject to an optional pitching moment constraint. 

• Effects of fuselage and nacelles be included in the design solution. 

• Constraints be applied to the design pressure distribution and/or 
local camber surface shape to provide physical realism. 

The design point options are discussed on page 49. 

The optimum loading selection is an extension of the methods used in 
references 9 and 11, through the addition of the configuration-dependent 
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loadings (due to fuselage and nacelles) and the addition of the pressure and 
ordinate constraint formulation. Also, more component loadings were added to 
provide flexibility in rearranging pressure distributions at the design point 
in order to meet the pressure constraint condition. The maximum number of 
loadings is 17, defined in the table on page 44. 

The wing design method actually consists of four solution steps: 

• Flat wing solution 

• Calculation of aerodynamic characteristics of selected component 
loadings. 

• Optimization of loading combination 

• Camber shape calculation 


Flat winq solution. - The first step in the wing design process is to obtain 
the flat wing solution for the given wing planform and Mach number, using the 
analysis form of the lifting pressure calculation described on page 53 . This 
solution is obtained to locate the wing aerodynamic center for use in 
computing the pitching moment at zero lift for the individual loadings; and 
also, because its drag-due-to-1 ift factor is an often-used reference point in 
judging twisted and cambered wing designs. 


Comp onent loadings characteristics . - For each selected component loading, 
section lift, drag, and pitching moment coefficients are evaluated at a series 
of span stations and integrated spanwise to obtain overall lift, drag, and 
pitching moment. In addition, the interference drag coefficients between 
loadings are computed, where the drag interference between loadings i and j is 
given by 



+ V AC (L* ,N* ) (4—0 (L* ,N* ) A(L»,N#) B(L*,N*) 

$S JUt p,j 3x i 

where the summations are carried out over the wing planform. For the 
component loadings defined over the planform by analytical expressions, the 
loading is scaled by the appropriate power of either semispan or average wing 
chord, so that the resultant component loading lift coefficient is 
approximately l.Q. 
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TABLE 1. 


DESCRIPTION OF WING LOADING TERMS 


Loading 

Number Definition 

1. Uniform 

2. Proportional to x, the distance from the leading edge 

3. Proportional to y, the distance from the wing centerline 

4. Proportional to 

5. Proportional to x^ 

6. Proportional to x(c - x), where c is local chord 

7. Proportional to x^ (1.5 c -x) 

8. Proportional to 2 (1 + 15 — F®*® 

c 

9. Proportional to ^2 ^ _ •] )2 

10. Elliptical spanwise, proportional to V (1 - W) 

11. Proportional to x, the distance from the leading edge of an 
arbitrarily defined region 

12. A camber-induced loading proportional to the body buoyancy 

loading 

13. A camber-induced loading proportional to the body upwash loading 

14. A camber-induced loading proportional to the nacelle buoyancy 

loading 

15. The body buoyancy loading 

16. The body upwash loading 

17. The nacelle buoyancy loading 
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Interference drag coefficients associated with the configuration- dependent 
loadings are special cases. The configuration-dependent loadings are those 
produced by fuselage upwash, fuselage volume asymmetry, and nacelles, as 
described on page 62. They are calculated at discrete points over the wing 
planform, and interpolated linearly for a complete definition where needed . 
Since there is no camber surface associated with the configuration-dependent 
loadings, the interference drag coefficient between the "regular" loadings and 
the configuration dependent loadings has only one term. To ease the 
computational task of evaluating component and interference characteristics, 
the configuration-dependent loadings are assigned a corresponding wing slope 
definition of zero. 

The nacelle configuration-dependent loading has yet another difference, since 
the flow fields of the "regular" loadings act on the fixed geometry of the 
nacelles, producing an appreciable axial force. The lifting pressure acting 
on the wing lower surface is assumed to propagate along Mach lines from the 
wing to the nacelle and the associated buoyancy force is computed. 

If Z constraints are used, the camber surface Z values of all loadings are 
interpolated at the planform locations where Z constraints are to be applied. 

These are used to constrain the optimum loading combination by requiring the 
specified Z value to occur at the constraint locations. 

If a fuselage is present, the wing solution proceeds outboard from the 
side-of-fuselage station, and carry-over lift and fuselage camber is used to 
define the force coefficients of the wing inboard of the side-of-fuselage. 


Fuselage contribution and "carry-over 11 lift . - If a fuselage is present in 

the solution, the lifting pressure distribution of the component loadings 
inboard of the side-of-fuselage station is of the "carry-over" type, rather 
than the analytically defined loading. I.e., the loading on the wing area 
occupied by the fuselage is a dependent function of the loading defined 
outboard of the side-of-fuselage station. 

The carry-over pressure distributions are calculated for each component 
loading (excluding the configuration - dependent loadings) using the analysis 
form of the lifting pressure solution described on page 54. The carry-over 
pressure distributions are applied to the fuselage camberline to obtain the 

corresponding drag terms. 

In addition, the fuselage contribution to configuration lift, drag, and 
pitching moment is included in the design point force coefficients. The 

fuselage coefficients are calculated in the lift analysis program for the 

input fuselage geometry as described on page 68 and transferred to the wing, 
design program. 
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Optimization of loading combination . - Given the component loadings and their 
interference drag terms, the aerodynamic characteristics of a wing supporting 
any combination of these component loadings can be calculated in terms of load 
strength factors assigned to each. The total lift coefficient for n 
loadings is 


i=n 


C L = 


i=l 


C L,i A i 


(29) 


where C|_ s -j denotes the loading and Aj its load strength factor. The total 
pitching moment coefficient at zero lift is C mo 


where 


C = 
mo 


i-n 


2 

•? =i 



(30) 


C s c 
mo,i M,i 

and the total drag coefficient is 


dC M c 
dC L L,i 


i«n 

J=n 

i-n 

2 

2 C D.i.J A i * 

T f C DWN,i A i 

i=l 

J=1 , 

i=l 


(31) 


where the terms ^ are the axial nacelle force coefficients for the 
component loadings. Tfie contribution of each loading to the lifting pressure 
coefficient, to the wing upper surface longitudinal pressure gradient, and to 
the camber ordinate at point j can be summed to give 


i=l 


(32) 


The load strength factors of the conf iguration-dependent loadings are 1.0, 
since those loadings are not variable. 
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Hence, from the component loadings data, the drag, lift, and pitching moment 
characteristics of an infinite number of wing designs can be immediately 
calculated in terms of the factors. The formal optimization of the wing 
loading involves the specification of the set of A-j's which gives least drag 
subject to the imposed constraints. 

According to the Lagrangian method, the solution for the optimum A-j values 
involves the system of linear equations shown in figure 4.4-4, where matrix 
notation is used for compactness. The term *C P) -j is the allowable lifting 
pressure coefficient corresponding to the ith pressure limiting constraint on 
the wing upper surface pressure. 

Partitions in the square matrix of the solution have been identified by number 
and letter in figure 4.4-4. The left, uppermost partition of this matrix 
(zone Al) and the adjacent row and column correspond to Grant's original 
design solution (reference 14). The next m rows (zones C1-C7) and columns 
(zones A3-G3) are introduced by m local constraints on the lifting pressure 
coefficient, which are added (if necessary) in applying the pressure level 
constraint, as described on page 50. The next q rows (zones D1-D7) and 
columns (zones A4-G4) are introduced by q constraints on the longitudinal 
gradient of wing upper surface pressure coefficient; these terms are added in 
applying the pressure gradient limiting condition. 

The next three rows (zones E1-E7) and columns (zones A5-G5) correspond to the 
constraints which set the loading factors A-j to 1.0 for the body buoyancy, 
body upwash, and nacelle buoyancy loadings. Each of these three rows and 
columns is present only if the corresponding configuration dependent loading 
is used. Zones F1-F7 and A6-F6 correspond to r constraints on ordinate, and 
appear only if ordinate constraints are actually specified by the user. The 
right-hand column (zones A7-G7) and bottom row (zones G1-G6) correspond to the 
constraint on pitching moment coefficient at zero lift; this constraint is 
also optional. 

The unknowns of the design solution are the n values of the loading factors 
Ai and the (m+q+r+5) Lagrange multipliers ( A-j)* These are contained in 
the left-hand-side column matrix. 

The seven partitions of the right-hand-side column matrix contain: either the 

negative sum of the nacelle interference drag produced by wing lift and 
fuselage carry-over drag, or zero if not used; the design lift coefficient 
less fuselage lift coefficient; the pressure level constraint values applied 
to the lifting pressure coefficient (these are generated automatically by the 
design computer program, if they are necessary); the pressure gradient 
constraint values (also applied automatically if they are necessary); the 
values of A-j for the three configuration-dependent loadings; the constrained 
values of camber-line ordinate; and the design pitching moment coefficient at 
zero lift less fuselage pitching moment at zero lift. 
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Notes: 1. Subscript co identifies carry-over 

lift contribution, 

2. Subscript f identifies 
fuselage contribution. 




FIGURE 4.4-4. -OPT/MIZA T/ON MA TRIX 



The maximum number of equations in figure 4.4-4 is 34. This number of 
equations can be reached if 17 loadings are combined, a pitching moment 
constraint is used, and the maximum number of constraints on lifting pressure 
coefficient are required. 


Design camber surface. - Given the set of factors, the design pressure 

distribution is known and the resulting camber surface is calculated using 
equation (27) on page 39. 


Desiqn point options. - Due to options in the wing design program to limit 
the solution extent (e.g., no C mo constraint or no ordinate constraints or 
no pressure constraints), some parts of the solution matrix on page 49 may not 
be used. However, for purposes of describing the solution options, it is 
assumed that all options are chosen. In that case, the wing loading solution 
proceeds as follows, using repeated applications of the corresponding matrix: 

(1) The wing loading having minimum lift-dependent drag with a. constraint 
on design lift coefficient is defined. This solution will include 
ordinate constraints if any have been specified. 

(2) A family of wing loadings having minimum lift-dependent drag with a 
constraint on design lift coefficient and a series of constraints on 
pitching moment coefficient at zero lift are defined. This series of 
solutions is presented in the form of drag due to lift factor, Kjr, 
versus C mo , and is referred to as the "bucket" plot. Ordinate 
constraints are not included. 

(3) The wing loading having minimum lift-dependent drag with constraints 

on design lift coefficient wing upper surface pressure and camberline 
ordinates is found. This solution may require a cyclic process for 
the pressure constraints that begins with (1) above. If (1) 
satisfies the pressure criteria on both gradient and level everywhere 
on the planform, then this loading is set equal to (1). If (1) does 
not satisfy the pressure limits everywhere, then a pressure 

constraint is applied at the wing planform location where the 
pressure gradient is most strongly violated, and a new solution for 
wing loading is found. Its pressure distribution is tested against 
the gradient limit, and another constraint is imposed if the limit is 
violated. Each cycle adds another pressure gradient constraint to 
the solution in addition to the lift coefficient constraint. When 
pressure gradient is satisfied, then pressure level is similarly 
tested against its limit. Constraints on pressure level are also 
added cyclically until solution pressure, level satisfies the level 
limit everywhere. Cycling continues until either the solution load 
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distribution everywhere satisfies the pressure gradient and level 
limit or until the number of solution constraints equals the maximum 
allowable number. In the latter case the limit level for pressure 
gradient is increased by 20%, and the cycle is started again. 

It has^ been found that pressure gradient and pressure level 
constraints imposed early in this cyclic solution can subsequently 
become unnecessary as other pressure constraints are imposed. 
Constraints unnecessary for satisfying pressure gradient and level 
criteria penalize the wing design unnecessarily by increasing drag. 
Such constraints can be detected by the algebraic signs of their 
solution Lagrange multipliers. Accordingly, a test for unnecessary 
pressure constraints has been added to the solution cycle. Such 
constraints are removed when detected and cycling continues. 

(4) The wing loading having minimum lift-dependent drag with constraints 
on design lift coefficient, camberline ordinates, and zero-lift 
pitching moment coefficient is defined. 

(5) The wing loading having minimum lift-dependent drag with constraints 
on design lift coefficient, pitching moment coefficient at zero lift, 
camberline ordinates, and wing upper surface pressure is found. The 
latter type of constraints are imposed, if necessary, in the cyclic 
fashion of (3). 

Examples of loading solutions (1), (2), (4), and (5) are shown in figure 

4.4-5. For this case, loading solution (I) satisfies the wing upper surface 
pressure constraint, and is therefore identical to loading solution (3). 

The effect of the _ number of loadings on the "bucket 11 plot and on loading 
solution (5) is illustrated in figure 4.4-6. Increasing the number of 
loadings from 3 to 11 results in a substantial theoretical drag decrease; 
however, this decrease is reduced substantially when pressure gradient 
constraints are added. 

Experience with the wing design program has shown an occasional "endless loop" 
fault in the process of adding and deleting pressure constraints under the 
constraint cycling described in paragraph (3). I.e., a given constraint will 

be successively added and then deleted and then added again. The problem 

occurs due to an improper solution matrix formulation caused by imprecise 
interference drag calculations between loadings. Imprecision in the 
interference drag calculations is inherent in the panel type approach, and it 
is not possible to predict when it will cause an endless loop in the 
optimization. It is possible, however, to detect the presence of the 
mathematically improper matrix formulation through an eigenvalue check and 
thereby avoid entering an endless loop. 

The test for a probable endless loop situation is to calculate the number of 
zero and negative eigenvalues of the solution matrix. The matrix solution is 
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Wing upper surface C constraint: Cp ^0.7 C p 


vacuum 



^mo 

FIGURE 4.4-5. — TYPES OF LOADING SOLUTIONS 
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NOTE: 

Flagged symbols denote 
solutions satisfying 



Maximum for 3 loading solution, ^ 

including thickness pressures. mo 

FIGURE 4.4-6.-WING DESIGN STUDY: CONSTRAINTS ON BOTH C mn AND UPPER 
SURFACE PRESSURE GRADIENT 
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unstable (probable endless loop) if there are any negative eigenvalues, or if 
the number of zero eigenvalues exceeds the number of imposed constraints. The 
number of imposed constraints is equal to (if used) one each for total lift 
and pitching moment, plus the number of fixed configuration-dependent 
loadings, plus the number of Z constraints, plus the number of pressure 
coefficient level and gradient constraints. If the eigenvalue test proves 
unstable, the program halts and an error message is printed.. The optimization 
requested must then be changed to permit a solution. (A discussion of such 
changes is given in the User's Manual, Part 2.) 


Lift Analysis Case 


Loading for a given camber surface . - Since the R function for the field 
point element is zero (the element generates no downwash or upwash upon 
itself), equation (27) can be rearranged to solve for the lifting pressure 
coefficient in terms of the field point slope and upstream influences: 


AC (L*,N») 
P 


i, dz 

4 C 

B dx 


(L*,N») 


(33) 


+ 2j*R(L* - L,N» -N) A(L,N) B(L,N) C(L,N) AC p (L,N) 


Therefore, the lifting pressure distribution, ACp, can be determined for a 
wing of arbitrary surface shape provided the calculations are performed in the 
proper sequence. The order of calculation of AC p (L*,N*) is from apex aft, 
and from the centerline to the right-hand wing tip for each L* row. In this 
fashion, all pressure coefficients within the Mach forecone from any element 
will have been previously obtained and no unknown pressure coefficients arise 
in the sunmation. The element fractions A, B, and C are as defined for 
equation (27). 


Theoretically ACn(L*, N*) defined by equation (33) is the pressure 

coefficient at the aft mid-point of the L*, N* element. The average pressure 
coefficient for the element, needed in subsequent calculations, is calculated 
by one of two alternative methods: 

1) An approximate average pressure coefficient is interpolated from the 
aft mid-point AC p value and the AC p value of the element 

immediately ahead. (No interpolation is performed if L*, N* is a 
leading edge element). This is the method used in the original 
formulation of the computer program. However, it was found that 

oscillations in the AC p values occurred near the wing leading 
edge. So, a 9 point smoothing equation was applied after unsmoothed 
pressures were computed for the entire wing, which essentially 
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removed the oscillations and brought the smoothed pressure 
distribution into good agreement with established results from other 
theoretical methods. 

The smoothing equation is of the form 


0.2A{L*-U)AC p (L*-U) ♦ 0.1*A(L*-3)AC p (L«-3) ♦ 0.6 a(L*-2)AC (L*-2) 

♦ 0.8A(L*-l)AC ( L*-l ) ♦ A{L*)AC (L») ♦ 0.8 a(L» +l)AC (L* ♦ l) 

P P P 

_/ ^ ♦ 0 . 6A ( L*+2 ) ACp ( L*+2 ) + 0.l»A{L*+3)AC p (L*+3) + 0.2A(L*+4)AC p (L«MO 

AC p (L*) « 0 .2A(L«-U) + 0.l4A(L*-3) + 0.6 a(L*- 2) + 0.8A(L*-l) ♦ A(L* ) + 0.8A(L*-H) 

♦ 0.6a(L # + 2) 0.tA(L*+3) + 0.2A(L»+U) 



where the A values account for element fractions. A is zero if the 
corresponding element is ahead of the wing leading edge, or aft of 
trailing edge in the case of a subsonic trailing edge. The 

supersonic trailing edge solution is treated as a special case, and 
the trailing edge is extended four elements to provide pressure 
coefficient data to fill out the smoothing equation. 


2) The other average AC p method is an alternate technique, which uses 
an aft element sensing approach. This method involves solving for 
preliminary AC p results for a given L*, N* element and the element 
immediately aft, then following it up with a second calculation to 
refine the preliminary results. The procedure is detailed in 

reference 12, but is summarized below and in figure 4.4-7. 


a) Calculate preliminary AC p values for a given L* = constant 
row. Designate as AC p#a (L*,N*)* 

b) Calculate preliminary AC p values for the element row 
immediately aft = L*+l, using AC p>a values. Designate as 

Ac p,b(L\N*)- 

c) Finalize AC p values for the original L* = constant row from one 
of the two following equations. For leading edge elements: 


AC (L*,N*) 

V 2 


1 + a(l«,n«) 

1 + A(L*,N*) 


AC (L*,N*) 
P»a 


A(L»,N») 

[l + A(L*,N*) 


AC V (L*,N*) 
P 


(35) 
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7 
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As 


r 

' — Fairi 


Fairing of preliminary results 
Integration of preliminary AC 


AC, 


— Final AC p (L*) 


FIGURE 4.4-7.-APPLICATION OF AFT-ELEMENT SMOOTHING TECHNIQUE (FROM REF. 12) 



For all other elements: 


AC P = | AC p a (L*,N*) + i. AC p;b (L*,N») 

The aft element sensing technique produces a substantial reduction in pressure 
coefficient oscillation. In cambered wing test cases run with the method, 
however, some oscillation was found to persist. Therefore, a 3 term smoothing 
equation was added (after all wing pressures are calculated), of the form: 


ACUL*) = °- 5 A(L*-1)AC p (L*- 1) + AC p (L*) + 0.5A(L*+1 )AC D (L*+1) (37) 

P 05A(L*-1) +1.0 + 0.5A(L*+1) 


In the computer program one of the two alternate smoothing methods is 
selected by an input code. Both produce essentially the same answers . When 
the pressure limiting option is used (discussed later), the aft-sensing 
smoothing technique is automatically selected. 

For the analysis solution, the pressure coefficients for all elements must be 
calculated. The force coefficients are, therefore, calculated from direct 
summations of loca'l pressures applied to each element, rather than employing a 
spanwise integration as in the design case. Lengthwise and spanwise lift 
distributions are obtained by summing the lift in the corresponding element 
rows. 

Leading edge suction force . - The drag summation may be obtained either with 
or without the leading edge suction force included (as requested by input 
codes). Both the maximum theoretical leading edge suction or the attainable 

leading edge suction force defined by Carlson in references 18 and 19 may be 

calculated. The referenced reports contain detail description of the methods 
involved, which are summarized here. 

The magnitude of the theoretical leading edge thrust developed by a thin 

lifting wing is dependent upon the upwash in the vicinity of the leading 
edge. This upwash determines the character of the associated pressure 
distribution. For a flat lifting wing with a subsonic leading edge, the 

pressure distribution near the leading edge is peaky, with the lifting 
pressure at the leading edge theoretically approaching infinity. The 
theoretical leading edge thrust is related to the limiting value at the 
leading edge of the parameter 


lim (AC p V?" ) = (ACAS 7- ) (38) 

x^O P 0 
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and the thrust per unit distance in the spanwise direction is 


t 


= -~-tan A ^1 - (3 2 cot 2 A (AC p 



( 39 ) 


The total theoretical thrust for the wing becomes 


C T 



c t dy 



b/2 

o 


7 r b 

Ts tan A 


yl- cot^ A (ACp 



(40) 


The central problem in calculating the leading edge thrust involves the 
evaluation of ( AC p • It is complicated by the fact (typical of 

panel aerodynamic methods) that the calculated pressure distributions are 
appprox imately correct but vary from "exact theory" values and are most 
inaccurate near the leading edge. In reference 18, Carlson performed a detail 
analysis of the nature of the errors associated with the presssure 
distributions near the leading edge for the. "Mach-box" numerical method. An 
adjustment parameter was developed (a shift in the x‘ location of the 
calculated pressure near the leading edge as a function of wing sweepback and 
Mach number) which significantly improved the pressure distribution details 
near the leading edge, with attendant improvement in the value of (AC p vx' j Q . 

The solution for ( AC p v^) Q is obtained by a least squares fit of AC p v 'x T 
versus x', using the fitting function shown below: 



Since it is known that the flat wing pressure distribution near the leading 
edge displays a square root singularity, the function 


ACpV^x 7 " - kj + k^x' 


(41) 


is used for the flat wing solution. Wings input as a cambered surface are 


57 



analyzed with the more general expression 


ACpVx 7- - k { 



( 42 ) 


All wings, flat or cambered, could have been treated with the cambered-wing 
fitting function, but flat wings were considered to be a special case because 
of some improvement in results and relatively minor programming complications. 

Leading edge sectional thrust values are computed for all spanwise rows of the 
wing grid system for a series of wing angles of attack. 

The sectional thrust values are resolved into drag or lift components 
according to assumptions of fully attached flow, fully separated (vortex) 
flow, or the attainable thrust definition of reference 19. The components 
associated with the fully attached or separated (Polhamus analogy) flow 
assumptions are illustrated in figure 4.4-8 for a typical leading edge segment 
having a local angle of attack of a- 6* The sectional data are integrated 
spanwise to get the total wing data for output in the force coefficient 
summaries. 


Attainable leading edge suction . - Carlson, in reference 19, further 

developed the leading edge suction solution through the concept of an 
attainable leading edge thrust. Based on a study of airfoil data and 

utilizing simple sweep theory, reference 19 derives an empirical estimate of 
attached flow and vortex flow components. The principal factors of importance 
were found to be wing leading edge normal Mach number and wing airfoil normal 
section parameters (leading edge radius, maximum thickness and position of 
maximum thickness). From these parameters and sectional Reynolds number, a 
factor K T relating the theoretical sectional thrust coefficient to the 

attainable attached-flow-type thrust coefficient is computed. The difference 
between the attached flow thrust coefficient and the theoretical sectional 
thrust coefficient is then treated as a vortex lift component, with associated 
drag and lift coefficient components. The attainable thrust force coefficient 
components for a typical wing leading edge segment are illustrated in 
figure 4.4-9. The attainable lift and drag sectional coefficients are also 
integrated spanwise and printed in the wing force coefficient summaries. 

It should be noted that the theoretical leading edge suction force is an 
approximate calculation due to the panel representation of the wing and 
corresponding pressure coefficient dependency. Particularly in the case of a 
cambered wing, the calculated leading edge thrust forces typically involve a 
small but significant inaccuracy and should be used with appreciation for 
possible error. 

Typical drag polars with the different leading edge suction options are 
illustrated in figure 4.4-10. 
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• Attached flow • Separated (vortex) flow; 

Polhamus analogy 




eg = c t sin (a - 5) 
c d = “ c t cos ^ 


c _ cos (a - 5) sign (C ) 

* M x = 0 


"d ' 


cos A 
c t 

cos A 


sin (o - 8) sign (C ) 
= 0 


FIGURE 4.4-8.— THEORETICAL LEADING-EDGE SUCTION LIFT AND DRAG COMPONENTS 
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• Attainable suction 
sensitive to— 

• Airfoil geometry 

• Norma! Mach number 

• Reynolds number 



Attainable thrust factor, K-j- = 


= 


sin (cos~1 Ky -5) 


cos A 


-1 


c g = -c t cos (cos K-j — 5 ) 



If K-j- = 1.0, attached flow 
Kj = 0., Polhamus analogy 


FIGURE 4.4-9.- ATTAINABLE LEADING-EDGE SUCTION ( CARLSON ) 
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Analysis of complete configuration . - The preceding discussion describes the 
theory used in the lifting pressure calculations for a wing camber surface at 
a selected angle of attack. The analysis program may also be used to 
calculate lifting pressure distributions and force coefficients for complete 
configurations over a range of angles of attack, adding in the effects of 
fuselage, nacelle, canard, and/or horizontal tail, as applicable. 

The program actually carries two solutions along: one for the. configuration 

at its input angle of attack, the other the incremental solution per degree 
angle of attack (called the flat wing solution). The interference terms 
associated with the two solutions acting on the other surface shape are also 
calculated. The summation of these effects into the drag polar and other force 
coefficients is performed by superposition, as described on page 70. 

Calculation of the complete configuration lifting pressure solution involves 
up to 7 principal tasks: 

• Isolated fuselage upwash field 

• Nacelle pressure field acting on wing 

t Pressure field due to asymmetrical fuselage volume 

• Wing/canard solution in presence of fuselage upwash field. Effects 

of the wing pressure field acting on nacelles are also calculated 

• Fuselage lift distribution in presence of wing downwash field 

• Horizontal tail solution in presence of fuselage and wing flow fields. 

• Superposition of solutions 

The calculation logic of these tasks is described in the following sections: 


Isolated fuselage upwash field . - Two alternate solutions for the fuselage 
upwash field are provided, selected by an input code. One is a slender body 
solution which is valid for a pointed body of circular cross-section shape. 

The second is a general slender body solution of arbitrary cross-section. 


Fuselage of circular cross-section . - The fuselage upwash (or downwash) 
field in the plane of the wing, canard, or horizontal tail is obtained from 
slender body theory. The fuselage must have a pointed nose and cannot have 
discontinuities in cross sectional area (e.g., a step in longitudinal area 
distributions). 

From reference 15, the velocity potential of a slender body of revolution, 
shown in figure 4.4-11, is: 


-g0(x, r,6) 


= 1 d A (x)n 
2ir dx 


x 


J_ a_ f 
2 2rr 3 x J 
o 


dA(x) 

dx 


iIn(x-X|)dxj + [a+e(x)] R^(x) 


cos d 
r 


(43) 
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Upwash propagates along Mach lines. 
Upwash at arbitrary point P is: 



U 


w 

u 


- JL l_ _ e (X) R 2 !*) h 2 -v 2 'l 

2ir dx ,2 r 2 V r 2 ) 

- -s!m 


Where:: 


u Q = upwash at input fuselage incidence (ct = 0) 

u = upwash increment per unit O’ 

A(x), R(x) = Fuselage cross-sectional area, radius 

6 (x) = Fuselage incidence due to camber 
(relative to reference plane) 

z,y,r - location of point at which 
upwash is calculated 

ol = angle of attack of reference plane 
U ~ free stream velocity 


FIGURE 4.4 - 1 1. -FUSELAGE UPWASH CALCULATION (CIRCULAR CROSS SECTION ) 
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The nondimensional radial velocity component is 


~ 0 r - — — - dA(x) - CL+ ...f( x ) R 2 (x) cos 6 

U 1 2ttt dx r 2 

and the tangential component is 

~Ur 4>e = -[a + e(x)l sin d 

r 2 

The vertical velocity component (upwash angle) is 

Tf = Tf 0 r cos 9 - <t>9 sin 8 


(44) 


(45) 


(46) 


Substituting for <£ r and and converting to Cartesian coordinates gives 

_g=JLdA£x)_L.- [a + 6(x) , R~(x) (z 2 -A < 47 > 

u 2lT dx f 2 r 2 \ r 2 / 

When the fuselage is at input incidence, the fuselage upwash is 


= _L dA(x) z _ e / x) R 2 (x) (z 2 - y 2 ^ (48) 

2tt dx r 2 r 2 \ r 2 ) 


and the incremental upwash per degree angle of attack increase is 


* = d w/u _ _ R 2 ( x) 
11 da r 2 



(49) 


Slender body theory assumes that the upwash field propagates normal to the 
body centerline. In the program, options are provided for the upwash field to 
be propagated either normal to the body centerline or along Mach lines (i.e., 
according to the characteristics of supersonic flow). 


i 
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Fuselage of arbitrary cross-section . - In cases where the fuselage 

cross-section is significantly different from circular, or where the fuselage 
nose is open, or where discontinuities in longitudinal area distribution 
occur, an alternative fuselage solution may be selected. The cross-sectional 
shape for this option is input as sets of Y and Z coordinates at a series of 
longitudinal (X) stations. 

The method used for the arbitrary cross-section fuselage solution has been 

referred to as general slender body theory. The computer program employed was 
developed by Bonner at Rockwell-International and was adapted for use in the 
present application. A description of the analysis is given in the Appendix 
and summarized here. 

The perturbation velocity components due to the fuselage are obtained by a 

linear combination of two solutions: 

• An axisymmetric solution (non-lifting) due to the thickness 

distribution of the equivalent body of revolution. The numerical 

method employs linearly varying line sources . and sinks. 
Discontinuities in the area distribution (defined in the input 
definition) are treated by step perturbation functions. 

• A cross-flow solution to account for fuselage angle of attack and 
cross-sectional shape. The method of conformal mapping of the known 
solutions for a circular contour is used for this analysis. 
Discontinuities in the fuselage longitudinal area distribution are 
treated by step perturbation functions. 

The cross-flow solution is obtained through the conformal transformation 
solution of Theodorsen and Garrick (reference 20). The method is quite 
general, requiring principally that the cross-section contour must be simply 
connected and describable by a single valued function in polar coordinates 
about the centroid of the section. 

Static pressures on the surface of the fuselage are computed and integrated 

over the surface contours to get isolated fuselage lift, drag, and pitching 
moment. Upwash flow angles in the plane of the wing are computed for the 
fuselage input angle of attack and per degree additional angle of attack and 
used in the wing/canard/tail solutions, as described on page 66. 


Nacelle pressure field acting on wing. - The pressure field caused by the 

nacelles on the wing is calculated by the Whitham solution, as described in 
Section 4.3. The nacelles are assumed to be bodies of revolution, and the 

pressure signature due to each nacelle is calculated for a series of spanwise 
stations from wing root to tip. The composite pressure signature at a given 

spanwise station is the sum of the individual nacelle signatures, with all 

pressure coefficients doubled to account for reflection from the wing. Either 
"wrap" or "glance" nacelles may be used, as discussed in Section 4.3. 
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Pressure field due to asymmetric fuselage volume . - Another term in the 
superposition approach is associated with the growth in fuselage area above 
and below . the wing. If the area growth is asyrrmetric (e.g., a low wing 
configuration), a differential pressure across the wing plane is created, 
which is additive to the effect created by the fuselage upwash field. 

An approximate method of computing the asymmetric fuselage pressure field is 
employed, using the area distributions shown in figure 4.4-12. To the actual 
forebody area distribution is added the growth in fuselage area in the wing 
region relative to the area at the leading edge of the wing-fuselage 
intersection. This is done for both the below- and above-wing area 
distributions. Pressure signatures for both body representations are computed 
using the Whitham technique at a series of spanwise stations, doubled to 
account for reflection from the wing, and differenced to get the lifting 
pressure distribution due to asymmetric fuselage volume. 

If the wing trailing edge is subsonic, the pressure fields above or below the 
wing revert to that corresponding to the actual fuselage area distribution in 
the region aft of a Mach line from the trailing edge of the wing-fuselage 
intersection. 

Use of the asymmetric fuselage calculation is controlled by an input code in 
the computer program. If the asymmetric calculation is not requested, the 
fuselage pressure field corresponding to a mid-wing arrangement is calculated, 
so that there will always be a thickness pressure field due to the fuselage, 
if present, for use in limiting pressure calculations (discussed later). 


Wing/Canard Solution in Presence of Fuselage and Nacelles 


The wing or canard . lifting pressure solution in the presence of fuselage 
and/or nacel les is performed as described for the isolated wing, with the 
following alterations: 

• The local surface angles of attack are increased by the fuselage 
upwash values for the purposes of computing pressure coefficients. 
(To compute drag, the pressure coefficient is applied to the local 
surface slopes, only). 

• In the region of the wing covered by the fuselage, the wing slopes 
are zeroed. _ This reflects the fact that the wing lift distribution 
in this region is of the "carry-over" type, only. (To compute drag, 
the carry- over lift is applied to the fuselage camber line shape.) 

• The canard pressure distributions are computed in the same fashion as 
the wing pressure distribution, with the canard also creating a 
downwash (or upwash) field for the wing. 
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FIGURE 4.4-12 FUSELAGE AREA REPRESENTATION FOR ASYMMETRIC 

FUSELAGE BUOYANCY FIELD 



The nacelle or asymmetric fuselage volume pressure fields are super imDosed 
upon the wing elements to obtain the effect of the nacelles or fuselage 
asymmetry on the wing,. If nacelles are present, the buoyancy drag of the wing 
lifting pressures acting on the nacelle area distribution is also calculated. 
This term is obtained by transferring local wing pressures aft along Mach 
lines to act upon the nacelle area growth, in a fashion paralleling the 
wi ng-th ickness-on-nacel le term of the near-field wave drag program. 

Fuselage force coefficients in presence of wing. - Two different fuselage 
force coefficient solutions may be obtained in the program, corresponding to 
the circular cross-section or arbitrary cross-section solutions described 
previously. 

IjJfselacje of circular cross-section . - From slender body theory, the fuselage 
lift distribution is given by 



c R(x) 


dA(x) a 

[a + e(x)] + A(x) — [e(x)] 

dx 


dx 


(50) 


using the notation of figure 4.4-11. 

The first term of the equation is the slender body term associated with a 
straight body at angle of attack; the second term is due to the curvature 
change associated with body camber. 

In computing the. fuselage .1 ift distribution, the downwash (or upwash) from the 
wing/canard is included in the local angle of attack term, e(x). For drag 
computations, the lift is applied only to the physical fuselage attitude. 

Fuselage force coefficients (lift, drag, and pitching moment) are calculated 
from the slender. body lift distribution and converted to a wing reference area 
basis for summation into the complete configuration characteristics. 


Fuselage of arbitrary cross-section . - Static pressures on the surface of the 
fuselage are computed and integrated over the surface contours to get fuselage 
lift, drag, and pitching moment. In computing the surface presssures, the 
downwash (or upwash) from the wing/canard is included in the local anqle of 
attack term. 

For the circular cross-section fuselage, the fuselage "camber drag" is 
computed by integrating the fuselage lift along the fuselage mean line, as 
noted above. For the arbitrary cross-section fuselage, there is no such 
meaningful Z reference line. For consistency in drag bookkeeping (and a 
direct comparison in case the cross-section is circular), the fuselage lift 
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along the Z centroid line is computed and printed. In addition, the 
integrated fuselage lift, drag, and pitching moment from the surface pressures 
are calculated for a series of fuselage angles of attack. The drag in this 
case consists of both drag due to thickness and drag due to angle of attack. 
An example of the fuselage data printout is given in Part 4 (Test Cases). 


Horizontal tail solution . - The horizontal tail lifting pressure distribution 
Ts calculated Tn the presence of the fuselage upwash field and wing/canard 
downwash field. The solution employs the same logic of partitioning the tail 
into exposed and carry-over regions used in the wing and canard pressure 
computations. 

Since the horizontal tail is assumed to be aft of the wing, the fuselage and 
wing/canard upwash (or downwash) field is first calcul ated, then the tail lift 
distribution and force coefficients are computed for the desired tail angles 
of attack. These are summed into the wing/fuselage, etc., solution for each 
tail angle, resulting in a set of force coefficients for each tail setting in 
the same fashion as wind tunnel data are obtained. 

Calculation of the wing/canard downwash field is performed by extending the 
wing grid system aft to include the horizontal tail. With all pressure 
coefficients on the wing or canard previously established, the effect of the 
wing/canard on the horizontal tail is obtained from 

Wing 

2 R (L* - L, N* - N) A (L, N) B(L, N) C(L, N) AC„ (L, N) 
r p (51) 

R (L* - L, N* - N) A(L, N) B(L, N) C(L, N) AC p (L, N) 


A W,C " 
Canard 

+ 

2 

T 



in which the wing, canard, and tail are all assumed to be located in 
essentially the same plane. 

Calculations of the horizontal tail lifting pressures are then performed using 
the analysis form of the lifting pressure equation 

AC (L*, N*)= & (L*, N*) 

F p dx 

1 Tail 

T 2 R (L* - L, N* - N ) A(L, N) B(L, N) C(L, N) AC p (L, N) + A w c 

with the wing/canard term added to the summation of upstream effects, 
the wing/canard pressure solution, the fuselage upwash is added 


(52) 


As in 
to the 
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physical tail slope for the purposes of computing lifting pressures, but 
removed for the calculation of the drag term. The carry-over lift of the 
horizontal tail is applied to the fuselage camber line to obtain the drag 
interference of the tail on the fuselage. 

The effect of wing downwash on the horizontal tail lifting pressures is 
significant. The data of figure 4.4-13 show the theoretical drag and pitching 
moment increments at constant total lift for a typical configuration, with and 
without wing downwash included. {Comparisons of the theoretical calculations 
with experimental data are presented in the User's Manual, Part 2.) 


Canard or horizontal tail downwash shift . - The lift analysis program contains 
a feature to shift the downwash field of wing or canard in the horizontal 
plane. The basic theoretical solution assumes the downwash field moves 
directly aft, instead of following the fuselage contour, as illustrated in 
figure 4.4-14. The shift feature can have significant influence on the 
calculated interference characteristics, as shown in figure 4.4-15. Program 
options allow for: 

• No shift 



• Shift according to the side-of-fuselage span stations of wing, 
canard, or tail. 

• Shift according to an input dimension. 


Superposition of Solutions 


The arrangement of terms in the superposition of the complete configuration 
force coefficients is: 


C L=2C La = () +a2 AC L 

(53) 

mm 

C m = 2 C m a -q i ~ aS Ac m 

(54) 

uJ 

]> 

& 

C D^C Da = o +K , (C L -C La = o ) + K 2 (C L -C La = o )2 

(55) 

J ' 

m 

3 

13 

* 

where the lift, drag, and pitching moment coefficients at a= 

0 correspond to 
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M = 2.7, C L = 0.10 


Coefficients based on wing area 
Increments, from tart on/tail off data 






Flow angle at tail 
downwash plane, deg 



FLOW ANGLE AT 
TYPICAL DOWNWASH PLANE 

Wing side of 
fuselage 


Upwash station 



station 


FIGURE 4.4- 14.- DOWNWASH PA TH CORRECTION 
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the individual fuselage, wing, tail, etc., solutions at input configuration 
incidence, and the AC[_ and AC m coefficients correspond to the 
incremental (flat wing) solutions per degree angle of attack. 

The drag coefficient equation includes the interference terms between the 
solution at input incidence and the incremental solution per degree. With 
subscripts C and F denoting input incidence and the flat wing terms, 
respecti vely, the drag coefficient equation is: 


C D “ C D C + ( C 


D FOC + Cd 


COF 


a 

a F 


C, 


(56) 


( C D, 


+ C 




or Cn , Cn l + -^o c. . ~ d cq f ) (Cl _ Cl , + AAALA. -V : 

C ^Lp 


D ” l D 


C 


T 


Li 


where a= configuration angle of attack 


a 


F= flat wing incidence, .01745 radian 


(57) 


C L C S 2 A C P C a E 


1 


: D r ~ S 2 A C P C a E a C 


C L F “ S 2 A C P P a E 

r F n u p r - Lp 


C Dr “ a F C Li 


with A £ = local element area corresponding to AC p 


S = reference area 


a. 


* 

c = local element slope 



C 


D 


FOC 


1 

— c 2 A Cd Ar CL 


E C 


interference drag of flat wing 
pressures acting on cambered 
wing slopes, per degree. 


C D = ap Ct interference drag of cambered 

L wing pressures acting on flat 

wing incidence, per degree. 


with the summations carried out over the wing planform. 

Therefore, for each of the solutions, involved in the superposition, it is 
necessary to calculate both an input incidence solution, a flat wing solution, 
and the interference drag terms between the two solutions. All of the force 
coefficients are referred to an input area, moment center, and moment 
reference length in the sunmations. 
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Limiting pressure solution . - Linear theory imposes no restrictions on the 
allowable wing pressure coefficients. Particularly in the case of subsonic 
leading edge wings, large upper surface pressure coefficients may be computed 
near the outboard wing leading edge that exceed (are more negative than) 
pressure coefficients corresponding to a vacuum. These pressure coefficients 
can occur at moderate wing angles of attack, due to the strength of the upwash 
generated by the inboard wing. 

Several investigators have examined experimental data to define the minimum 
attainable upper surface pressure coefficients. One such correlation is shown 
in figure 4.4-16 (from reference 16), which indicates a physical pressure 
coefficient limit of approximately .8 vacuum. The mechanism of this limiting 
is associated with local leading edge flow separation. However, the 
phenomenon may be approximated in the wing analysis program by limiting the 
calculated linear theory pressure coefficients to a specified fraction of 
vacuum Cp. 

In the analysis program, pressure limiting is an optional feature, controlled 
by an input code. If limiting is requested, then a set of configuration 
angles of attack is required for the solution, since superposition will not 
apply after limiting occurs. Also, a definition of the wing thickness 
pressures is required (transferred over from the near field wave drag program 
module), since it is the sum of lifting plus thickness pressures that is 
limited, rather than the lifting pressure alone. 

Pressure coefficients, as calculated, are separated into upper and lower 
surface values. Thickness pressures (wing plus fuselage) are added and the 
upper surface pressure coefficient tested against the limit pressure 
coefficient. If the limit is exceeded, the calculated lifting pressure 
coefficient is adjusted to a level which causes the total upper surface 
pressure coefficient to match the limit value. 

Comparisons of the "limiter" calculations with measured force and pressure 
coefficient data are presented in figures 4.4-17 through 4.4-21 for 2 arrow 
wings of reference 1. Both were 70° sweptback wings, having 3 percent 
biconvex airfoils, and wing design lift coefficients of 0 (flat) and .08. 

The basic thickness pressure comparison is shown in figure 4.4-17. 
Comparisons of upper and lower surface pressure coefficients at lifting 
conditions of C L = .13 and .25 are shown in figures 4.4-18 and 4.4-19. In 
the case of the flat wing (figure 4.4-18), several fractions of vacuum 
pressure coefficient are illustrated; the design C|_ = .08 wing data are 
shown for .7 vacuum limiting only. 

The limiter feature exhibits considerable improvement in the detail pressure 
coefficient comparisons over unlimited linear theory as Cj_ increases. 

Force coefficient comparisons for the two wings are shown in figures 4.4-20 
and 4.4-21. Pressure limiting at approximately .7 vacuum improves the 

comparisons appreciably at the higher C|_ values. 
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FIGURE 4.4- 17.— THICKNESS PRESSURE COMPARISON 
(C p DUE TO THICKNESS ONL Y) 
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M = 2.05 


Numerical method 





Experiment 

(NASA TN D-1264, wing 1) 

O Upper surface 
□ Lower surface 


V/"2 = 0.5 




x/c x/c x/c 


No C p limit 

Limit at 0.7 vacuum (C p = -0.24) 

Limit at 0.5 vacuum (C p = -0.17) 

Limit at vacuum (Cp = -0.34) 



C, =0.127 
L exp 



x/c 


x/c 


FIGURE 4.4-18 -PRESSURE COEFFICIENT COMPARISON, WING 1 (FLAT WING) 









M = 2.05 




Experiment 


/ 

(NASA TN D-1264, wing 2) 

Numerical method 


O Upper surface 

No C p limit 

— ' / 

□ Lower surface 

Limit at 0.7 vacuum (C 


.2 
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0 .2 .4 .6 .8 1.0 0 .2 .4 .6 .8 1.0 0 .2 .4 .6 .8 1.0 0 .2 .4 .6 .8 1.0 0 .2 .4 .6 .8 1.0 

x/c x/c x/c x/c x/c 

FIGURE 4.4-19.-PRESSURE COEFFICIENT COMPARISON, WING 2 

(TWISTED AND CAMBERED WING, DESIGN C, = .08) 















5.0 INTERACTIVE GRAPHICS 



Interactive graphics for use with the design and analysis system are optional, 
employing the NASA-LRC cathode ray tube (CRT) display and associated 
software. Use of the graphics option is requested by an executive control 
card (described in User's Manual, Part 2). 

The principal uses of the graphics routines in the design and analysis system 
are to display the configuration, edit input geometry, and to display and/or 
alter the basic program calculations. Limited capability to redirect the 
system calculation sequence is available from the CRT console during program 
execution. 

The CRT display and program coding for the interactive graphics setup are 
based on the NASA-LRC system. However, all display portions of the coding are 
subroutined or overlaid from the basic programs, so that the system could be 
readily converted to other CRT arrangements. 

Details of the interactive graphics portion of the design and analysis system 
are given in (User's Manual, Part 2). 
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APPENDIX 



General Slender Body Theory 


The solution for a fuselage of arbitrary cross-sectional shape employs a 
computer program developed by E. Bonner at the Los Angeles Aircraft Division 
of Rockwell International. The computer program is described in detail in 
Rockwell report NA-74-687, which is excerpted here to supplement the theory 
description on page 65 of this report. Nomenclature in this section, which is 
taken from the Rockwell report, is presented on page 93. 

Theoretical approach . - The approach used has been referred to in the 
literature as general slender body theory. The solution is based on limiting 
(slender body) solutions to the three' dimensional non-linear potential small 
perturbation equations of motion. 

[ 1 -Mi-(r+ 1)M£ /U-] 0 XX + 0 yy + 4> zz = 0. 

fl) 

with bounding conditions v ' 

* 0 on S 

on 

The perturbation velocity potential on the surface and in the neighborhood of 
a slender three-dimensional system satisfying equation (1) is in the form 


0 ( x, y, z, Moc) = 0 2 ( y, z; x) + g ( x, M~) 
for , 

VMJ-l * <<i 


( 2 ) 


$2 is the solution of Laplace's equation 


e 

R C 


(Myy+(M ZZ =0 


(2a) 


satisfying three dimensional boundary conditions in the YZ plane at the axial 
station X and g is a contribution dependent on Mach number and the equivalent 
body of revolution but not on the shape of the contour. The crossflow 
solution ^ 2 1S in general comprised of the following three terms: 


02 ( y> Z, x) 02 f t 02 ,Cl+ 02 p 


(3) 


where ^2 t 1S ^ w0 di m ensional incompressible flow solution for expansion 
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or contraction of the section (thickness problem) and 4>'£ a and 4>2 » 3 are 
the solutions for angle of attack or sideslip, respectively.’ The method of 
conformal mapping of the known solutions for a circular contour is used for 
this analysis. 

The solutions for a non-lifting body of revolution are used to establish the 
function g appearing in equation 2: 


g ( x, Me) r x ™ q [ 0g ( x, r, M») - 02 ? B ( y, z; x) ] 


where 02, B = /2TT S' (x) In r and 4>g are the crossflow thickness and 
total thickness solution, respectively. The latter term is established from 
the more simple problem of compressible axisymmetric flow 



[ 1 - M2, - Me - 0 X ] 0 XX + 0 r / r + 0 rr = 0 


(5) 


The three dimensional perturbation potential is obtained from a superposition 
of the crossflow and axisymmetric solution. 


0 - 02 , t + 02 , /3 - 02 , B + 0B 


( 6 ) 


The body axes Cartesian velocity components u, v, w follow from the standard 
potential differentations 0 X , 0^, and <J) Z , respectively. 


Crossflow solution . - The approach here is to transform the known complex 
potential solutions (ref. 21, p. 29) due to thickness and attitude for a 
circular contour with the center at the origin to the boundary of interest 
using conformal mapping theory. (The transformed boundary encompasses the 
class of simply connected contours which are describable by a single valued ) ~ m 

function in polar coordinates about a centered origin.) -Vl 


W 2 b = rv r In a = S (x) In q 

’ 27T 


W 2 , ju = - iv c e_i ^ (° ~ R o e 2i/i to) 


(7) 


) S 

Di 

) 3 

Q I 


( 8 ) 



02 + i 02 = w [ a ( f) ] = W (O') 


(9) 
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C and o are the complex variables in the physical and transformed plane, 
respectively, as shown in the sketch below: 



oo 

f = 0 + a n (x)+^li^+ a -200 + ...» 2 -Pi*) 

a n=-l c n 


The evaluation of the transformation constants a n is accomplished according 
to the method of Theodorsen and Garrick (ref. 20) by casting the mapping into 
the equivalent form 


SC n /a n 
o e t n 11 


( 11 ) 



The constants C n are related to the physical contour by the following set of 
integral equations: 


pH ?2 ti 

C n = A n + i B n = ~J 0 e 1 na In R dtp 


( 12 ) 


87 




where 


1 r> 2?r 

In R 0 = f In R d 0 

J 0 (13) 

and 

$ ~ + Z [ sin n 0 - cos n 0] 

n= 1 l R? R? J 

(14) 

This set of expressions is solved by the method of successive approximations 
using the known function R(-O-) describing the actual contour as initial data. 

The original set of mapping constants a n of equation (10) are related to the 
constants C n by 



a 0 - Cj 

ai = C, ^/2 + C 2 


a n - 1 = — ,E, 

n n - 1 


a , . C , + C 
n-n -In n 


(15) 


Contour points in the physical^ .(Ri » ^j) plane are located in the 

transformed plane by a=Ro e ’ where 0-j (0-j ) is given by 

equation (14). Transformed field points are established by rearranging 
equation (10) in the following successive approximation form: 
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It has been established that the transformation procedure is satisfactory for 
a wider variety of "fuselage-1 ilce" contours using 25 transformation terms or 
less (n^-25). 

The crossflow solution perturbation velocity components due to thickness and 
angle of attack are 


^2,t_=JL R.P. [S\x) J_ + S'" (x) In a] 
Uco 27 r o dx 


v 2,t = S R.P. 

1 

da 

2tt 

a 

"df 

w 2£_ ~ 

.. S\x) 

I.P. 

1. da 

Uoo 

2tt 


a df 

u 2,a = I.P. 

^ 

o 

to 

+ 

j 

-2R 0 

UooU-Cx) 

l ° 


dx 

v 

2, a 

= IP. [l 

C*o) 

2 +l]d£ 


U^cRx) a d£ 


dR c dRo 1 
dx a 


(18) 


dZp + J_ diz Q 
dx oCfx) dx 2 a 

(19) 


W2 - a — = R.P. [ (^.o) 2 +l] -1 

Uooa(x) a df 

where a(x) = a-dz Q /d x 

Axisymmetric solution. - The linearized thickness solution is numerically 
evaluated by the superposition of K linearly varying finite line sources and 
sinks. 



4>B (x, r, 


x-j3 r 

= f TOd£ 

-60 ^(x-f) 2 -/3 2 r 2 


( 20 ) 
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where 


f(*) = 


k 

£ 


k= 1 


a k 


? k <l<2 


thus 


k 

0 B (x, r. M 00 ) = £ 
k = 1 



and 


x - Br 

<P k (x, r, MJ = / 
x k 


~ a k ^ ~ x k> d I 
y (x - £)“ - (3 2 r 2 


(23) 


The unit velocity components in cylindrical coordinates 


u k (x,4) - d(p k /dx 


v k (x,r) = 30 k /8r 


(24) 


are given by 


-1 


u k (x,r) = - cosh xj-x k ^-ln 


2 x - Xi 


P r 


k (x,r) = (3 


Itt 1 ) 


2Lz_ x k 

r 


(25) 


3 

. © 

1H 

sm 

3 
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The axisynmetric solution velocities are obtained from the sum of the 
individual line contributions 


U B 


(x, r) = K 

l 


1 


a k u k < x > r ) 


a4>B 

3x 


V B 


(x, r) = K 

l 

k = 1 


a k v k ( x > r ) = d$B 

dr 


( 26 ) 


The unknown line source constraints, a K , are established by satisfying the 
flow tangency condition at the surface of the body through solving the 
associated set of simultaneous equations. 


Complete solution . - The total perturbation velocity components are obtained 
by a linear combination of the crossflow and axisymmetric solutions. 


+ u 2,t» aW + u 2,ff g - s'(x) InR 

Uqo Ut» a(x) (SUqo 2n 


3$ 

dx 


-1 


+ 


Ug 

U 


(x, r) 


oo 


= V' 


—L + V " ,CL + v 2,/3 ^ - s'(x) cosd + Vg (x, r) cos'd 


U, 


U, 


a(x) U c 


£U C 


2ir 


R 


U r 


b(j) 

Ty 


— + /3 


( 27 ) 


w = w 2 t + w 2>C( + w 2 (3 P . S’(x) sine- + v B (x, r) sin e 

Uoo oWu*, 2 tt R U„ 

3 (p - ex. 

3z 

Where the last term on the right hand side is the equivalent body of 
revolution surface result. 
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The static pressure coefficient is established from these velocities using the 
body axis small perturbation form of Bernoulli's equation: 
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NOMENCLATURE 


a k 


f (0 


9 

I.P. 

M 


k th 

line source coefficient, equation (23) 

Conformal transformation complex constants, equations (10), (l5) 
Conformal transformation complex constants, equation (l2) 
Static pressure coefficient 
Source strength per unit length 
Equation (4) 

Imaginary part 
Freestream Mach number 


R 

R.P. 

R o 

S 

S' 

s 

u, v, w 


Radius of equivalent body 
Real part 

Contour radius in transformed plane, equation (13) and 
figure on page 87 

Equivalent body cross sectional area 

dS 

dx 

Arbitrary surface 

Body axis longitudinal, lateral, and vertical perturbation 
velocities 


U B> V B 

u k> v k 

u 2 

v 2 

w 2 

u 2,t’ v 2,t> w 2,t 


Equivalent body axial and radial perturbation velocities, 
equation (29) 

k th 

line source unit axial and radial perturbation 
velocity, equations (27), (28) 

f*2 = R.P. 

OX dx 

= R.P. ^ -- 

3y 3a d£ 

—2 _ T p dw da 

da d£ 

Crossflow perturbation velocities due to thickness (growth) 
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U2>a, V2,a, 
W2,a 

u 2 > 0 > v 2 > 3 , 
W2» 0 

Uco 

v c 

w 

w 2, B 

w 2>/- 
x, r, ■©• 

x, r-j , <J) 
x, y, z 

x k 

d£_ 

do 

do 

dx 

z o 

a 

a(x) 

6 

Y 

C 

M 

5 

a 


Cross flow perturbation velocities due to angle of attack 
(vertical translation) 

Crossflow perturbation velocities due to sideslip (lateral 
translation) 

Freestream velocity 

Cross flow velocity, f a 2 (x) + p 2 
Physical plane complex potential, 4 > 2 + i^2 
Transformed plane complex potential due to thickness 

Transformed plane complex potential due to attitude 

Physical plane cylindrical coordinate system 
Transformed plane cylindrical coordinate system 

Right-hand body axis Cartesian coordinate system 
k th 

line source origin 

1 °° a 77 

mapping derivative (equation (18)) = 1 - — 2 V — 

0 7 ]= 1 

mapping derivative (equation (18)) = - 2 — da? ? 

df 77=0 o ^ dx 

z location of conformal centroid, I.P. a Q 
Angle of attack 
a-dz 0 /dx 

/m 2 -1 or sideslip angle 
Ratio of specific heat 

Complex variable in physical plane, y + iz 
See figure on page 87 
Dumny axial variable 

Complex variable in transformed plane, y + iz^ 
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4 > 

4>b 

4>2 

4>2, B 
4>2, t 
4>2,ot 
4 > 2,8 

$ 

subscripts: 

B 

c 

k 

X, y, z, r 
xx, yy, zz, rr 


Total perturbation velocity potential* 

Axisynmetric perturbation potential 
Crossflow potential 

4>2,t f° r circular contour, equation (7) 

Crossflow perturbation potential due to thickness 
Crossflow perturbation potential due to angle of attack 
Crossflow perturbation potential due to sideslip 

Total velocity potential, Uoo [ x - /3y + az] + <p 

Limiting (slender body) result 

Equivalent body 
Centroid of contour 
k th 

line source 

a/a x , a/a y , a/a z? a/a r 

3-/dx~, d^/dy~, d^jdx 9^/9 


* All potentials are dimensional. 
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